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DRUG DISCOVERY USING MULTIPLE MEMBRANE MIMETIC AFFINITIES 

Field of the Invention 

This invention relates to piediction of biological properties More 
particularly, the present invention is directed to the measurement of multiple 
membrane affinities of test compounds, methods and compositions useful for acquiring 
data characteristic of such affinities, and a method and system for using such data 
alone or in combination with other molecular descriptors for the prediction of 
biological activity. 

Ba ckground and Summary of the Invention 

There has been much research and development effort directed to the 
definition of screening methods capable of identifying drug leads. The goal of such 
efforts is to define efficient methodologies for predicting biological activity in vivo by 
using empirically definable (or calculatable) descriptors to predict biological activity 
without the time investment and expense in costly animal studies. For many years 
most drug screening methods were based on conventional biological activity assays. 
More recently assays have been developed that predict cell membrane transport 
properties of test compounds Generally drug leads are generated by comparing 
biological and physical properties of test compounds with known compounds having 
recognized biological activity in vivo There is a significant body of literature directed 
to prediction of biological activities based on comparison of physical, chemical and 
biological descriptors and the use of pattern recognition analysis of such descriptors as 
part of drug screening protocols 



compounds for probable biological properties based principally on correlation of 
numerical values characteristic of their interaction with two or more membrane 
mimetic surfaces with corresponding values for control compounds of known 



The present invention is directed to a method of screening test 



biological activities/function The meth^r! ic em>:Mfv 



: r onen ips ^nri nr pip!0£!C2! aci;v;;;c^ ! he mem bi due binding properties ot tesi 
compounds can be calculated, or they can be determined empirically with use of, for 
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example, liposomes, immobilized artificial membranes, such as those described in U.S. 
Patent 4,93 1,498, the disclosure of which is incorporated herein by reference, 
Langmuir Blodget films, computer chips or similar devices with immobilized lipids, 
capillary zone electrophoresis columns coated with membrane lipids, and the like 
In one embodiment of this invention vector calculus is utilized to 
pattern match membrane interaction values of test compounds with control compounds 
to predict biological properties. The pattern matching protocol can be applied to data 
sets containing only values characteristic of multiple membrane interactions, or such 
data sets can include other biologically significant molecular descriptors such as 
molecular surface area, molecular weight, dipole moment, octanol-water partition 
coefficients, molecular volume, membrane diffusion coefficients, metabolism rates, cell 
efflux rates, etc. 

In another embodiment, membrane binding data are obtained for test 
compounds and control compounds for use in accordance with this invention using 
immobilized artificial membrane chromatographic substrates in high pressure liquid 
chromatographic systems using aqueous mobile phases. Data relevant to the 
thermodynamics and kinetics of compound/membrane interaction is reflected in 
retention time and peak width, respectively. All data are preferably normalized relative 
to a standard compound or a set of compounds, for example, a set of compounds 
having a common biological activity or function. 

This invention also provides novel carboxyl-functional, head group- 
protected phospholipids useful for preparing immobilized artificial membrane 
structures useful for acquiring membrane interaction data They are prepared by novel 
high yielding transphosphitidylation of phosphatidylcholine derivatives using 
phospholipase D in the presence of protected alcohols 

Brief Description of the Drawings 

Fig 1 is a 3-dimensional vector plot of membrane affinity data for 
compounds known to act at serotonin recer^^ 

■,ripouna> Miuwn to a^i a\ nopamine receptors 
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Detailed Description of th e Invention 

In accordance with this invention in vitro membrane binding model was 
established to predict biological activity of compounds independent of both receptor 
binding assays and in vivo experiments The membrane binding model was developed 
5 based on measurement of membrane binding constants of biologically active 

compounds on membrane mimetic surfaces. Collectively, the membrane binding data 
demonstrate that compounds within therapeutic classes exhibit similar binding profiles 
on amphiphilic/membrane mimetic surfaces and concomitantly similar tendencies to 
distribute into tissue lipid pools where activity is elicited. 

10 Lipid heterogeneity in cell membranes functions to regulate the 

distribution and activity of many drugs A drag compound with affinity for a particular 
lipid will accumulate in lipid pools enriched with that lipid. Examples include 
vinblastin and chlorpromazine which accumulate in the inner leaflet of plasma 
membranes because phosphatidyl serine (PS) and phosphatidyl inositol (PI) are 

15 enriched on the cytoplasmic side of cells Drug-lipid interactions have been used to 
explain biological activity, the volume of distribution of drugs, drug binding to 
membrane receptors, and drug conformation in membranes. Boundary lipids near 
membrane receptors can play a key role in sequestering compounds and presenting 
membrane associated compounds to the receptor, or they can participate directly in 

20 regulating receptor function. Although such may not be true for all membrane 
receptors, it has been documented for several receptors The results from the 
experimental work upon which this invention is based demonstrates that compounds 
with similar biological activity exhibit similar membrane binding properties. 

Due to chemically precise ligand-receptor binding criteria, the drug 

25 discovery process traditionally focused on the structure of compounds as the dominant 
factor governing biological activity. However, it has now been observed that 
structurally diverse compounds within a given therapeutic class exhibit similar 
membrane binding properties indicates a cooperative role between the receptor and its 
adjacent membrane Thu^ the m vitrr* ^. n ^ r ^^ n u; n .^: n ., ^ j 

impound tMicit a giver r-ioiogicai actrvii) omce oom ine receptor and its 
surrounding membranes contribute to a compound eliciting a particular biological 
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activity, receptor binding assays used in conjunction with the present in vitro 
membrane binding model provides a synergistic approach to screening potential 
therapeutic compounds. 

In accordance with one embodiment of the present invention, a method 
5 of screening test compounds for probable biological properties is provided. The 
method comprises the steps of identifying two or more membrane mimetic surfaces 
each having a unique amphiphilic composition. A set of control compounds is selected 
for comparison purposes. Each control compound has a known biological property, 
for example, a biological activity or interaction with a known receptor type. 

10 The control compounds can have a common property, or they can be 

selected to represent widely variant biological properties. For each control compound 
there is defined an ordered set of numerical values characterizing a biologically 
relevant interaction (e.g., affinity) of that compound with each of the selected 
membrane mimetic surfaces. The ordered set of numerical values for each control 

15 compound or each set of control compounds (i.e., a "training set" as described below) 
can be represented by the expression (C,, C 2 C n ) wherein n is the number of 
membrane mimetic surfaces identified and used in the screening method. A similar 
ordered set of numerical values <T l9 T 2 T n > for each test compound characteristic of 
its biologically relevant interaction with each of the respective membrane mimetic 

20 surfaces is determined The set of numerical values for the test compound is then 
compared with the sets of respective values for the control compounds, and the 
biological properties of those control compound having ordered sets of numerical 
values best matching the respective numerical values in the ordered set of values for 
the test compound is identified The values Cj, C 2 C n and T 1? T 2 T n can be 

25 determined by computer calculations or empirically by any one or more of a wide 
variety of art-recognized techniques for evaluating membrane interactions. 

In one embodiment of this invention the numerical values characteristic 
of membrane affinity are determined chromatographically using an aqueous mobile 
phase and a stationary phase comnrisinc a membrane mirne^"r ^-^-o c 

oaten; ^ ^ i express 1 ^ incorporated Herein p\ rcicrcn^t: 1 nc Lei in memo r an c 

mimetic surface'' as used in describing and defining the present invention refers to any 
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surface bearing immobilized amphophilic molecules (i.e., those having both lipophilic 
and hydrophilic portions capable of exhibiting some selective affinity for or otherwise 
interacting with a solute (e.g., a test or control compound) in a fluid phase in contact 
with the surface. The term is intended to encompass a broad scope of commercially 
5 available stationary phases detailed for use chromatographic applications. Preferred 
membrane mimetic surfaces are those described in the above-incorporated U.S. Patent 
No. 4,931,498. 

Thus, in application of one embodiment of this invention membrane 
binding properties of a test compound of unknown biological activity are compared to 

10 the membrane binding properties of compounds having known in vivo biological 
activity to assess the probability that the test compound will exhibit one or more 
biological activities in vivo The analysis is made by pattern matching membrane 
binding constants measured on immobilized artificial membrane columns in a high 
performance liquid chromatography system. Two or more immobilized artificial 

15 membrane chromatographic columns are used Each test compound is injected and 
chromatographed on each column and the peak width and the peak time of the test 
compound is recorded and normalized to an internal or external standard. Thus 
experimentally measured selectivity values (a k . = k'/k' standard) and normalized 
standard deviations (0^= o/ostandard) were calculated for both the test compounds 

20 and compounds of known biological properties Pattern matching using vector 

calculus, multivariate analysis or principal component analysis of the ak' and cc^ values 
of the test compound and the control compounds allows comparison of the membrane 
binding properties of the test compounds and each of the control compounds or, if the 
control compounds all have a common biological activity/property, average or mean 

25 membrane binding values of the set of control compounds for each membrane surface. 

Important criteria for optimum implementation of the use of 
immobilized artificial membrane substrates for providing reliable values characteristic 
of membrane interaction include use of a stable membrane mimetic surface, preferably 
the same or a similar mobile phase composition an- 1 c.i-n~ 1 ^ — j 

-dng^ : n measurement:* using trie mcmr-ranc sunacc ii vjjh ua> -io-ja\ an j aiso iui-iu- 
lot variation inherent in membrane mimetic preparations 
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The external standard should be stable, it should be soluble in the 
aqueous mobile phase (typically organic phase modified), and it should have a 
reasonable retention time using said mobile phase and a high UV extinction coefficient 
The external standard is typically injected in the HPLC/IAM system frequently to 
5 reflect temporal fluctuation in retention time. Normalization of the data eliminates 
variation among synthetic lots and from column deterioration and other factors that 
can cause temporal variability in retention times and peak width. One external 
standard which has been found to give good results when the mobile phase is 85% 
0.01M PBS (pH 7 4)/15% acetyl nitrile is 4-methyl anisole. In one embodiment one or 

10 more sets of compounds, the members of each set having a common biological activity 
or function are combined with one or more test compounds and ioaded as a mixture 
onto a liquid chromatograph-mass spectrometer system to provide a direct comparison 
of binding properties of test and control compounds in that chromatographic system 
Membrane binding constants (k 1 values) are calculated in accordance 

1 5 with the formula k' = (t r - t 0 )/t 0 wherein t, is the retention time and t 0 is the dead time 
obtained from each injection, however, t 0 in the denominator was obtained by 
averaging the experimentally determined t 0 from all injections for each day the data are 
collected. 

The use of chromatography on immobilized artificial membranes 
20 provides data characteristic of both equilibrium binding to the respective membrane 
mimetic surface (k* values) and the kinetics of mass transfer between the immobilized 
artificial membrane surface and the mobile phase (i.e. 6 values). Similar to all 
chromatography columns, columns utilizing immobilized artificial membranes exhibit 
column-to-column and lot-to-lot variation This unavoidable variation in the 
25 experimentally determined k' and 6 values for each test compound and each control 

compound is eliminated, for pattern matching purposes, by normalizing k* and 6 values 
to an internal or external standard compound chromatographed on each immobilized 
artificial membrane column Thus k' and 6 values for the test compounds are each 
divided bv the respective k' ^nH v^m OC ^r*h n ^^^^^ ♦ 11 ; 

i^n.urs TM-j s>j!^-t!v::\ , ;nc jmirioDih/eu auiiiciai membrane chromatographic surlace 
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for molecular recognition of the analyte, a is a thermodynamic property and 
corresponds to a ratio of membrane partition coefficients. 

The values T„ T 2 , T n (a k . or a^) for each test compound and the 
corresponding values C u C 2 , C n for each control compound can be conveniently 
stored in any database management system for subsequent retrieval and analysis. The 
comparison of the set of numerical values for each test compound with the sets of 
respective values or the control compounds can be carried out using any of a wide 
variety of mathematical analytical techniques. 

In one embodiment of this invention the comparison is accomplished by 
multidimensional vector analysis of the data wherein the ordered set of numerical 
values for each test compound is compared to those corresponding values fur each 
control compound, or, where the control compounds are all selected to have a 
common biological activity or biological function, an average or mean values for the 
selected set of control compounds. Thus, in one embodiment the control compounds 
best matching the values for the test compounds are those for which the angle 6 in the 
formula cosine 0 = (TjC, + T 2 C 2 + . . . T^J/OV + T 2 2 + . . .T n 2 )y 2 (C, 2 + C 2 2 . . 
.C n 2 y/i is less than about 20°. The best match of membrane binding patterns between 
the test compound and the control compounds is that where the angle 0 in the above 
formula is a minimum Where the number of membrane surfaces used in the method is 
two or three, each ordered set of numerical values for the test compound and the 
values for each respective control compound can be displayed graphically as a vector 
in a 2- or 3-dimensionaI coordinate system, respectively, to facilitate comparison of the 
numerical values of the test compound with those of the control compounds 

In one embodiment the test compound data are displayed as a uniquely 
discernible vector quantity, i.e., in a color different than the control compounds 
Where the control compounds comprise a set of compounds having a predefined 
biological property, the numerical values for each member of said set of control 
compounds can be displayed as vector quantities uniquely discernible as a member of 
said set of control compounds and visiblv distinguishable fr,^ *hn 

, t iL'Jeimt'j Dioiomcdi propem r.^n nr mathematically rnanipuldieu to define a mean 
or average membrane interaction value for each membrane species The mean or 
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average vector quantities calculated for the set of control compounds can be used in 
the above-described pattern matching analysis or they can be displayed graphically as 
uniquely discernible vector quantities with those of one or more test compounds 
Thus, in one embodiment of this invention the membrane binding 

5 properties of six hallucinogens were measured on immobilized artificial membranes 
(IAMs) prepared from phosphatidylcholine (IAM.PC), phosphatidylserine (IAM.PS), 
phosphatidylethanolamine (IAM.PE), and sphingomyelin (IAM.SM) Mean vectors 
were calculated for several therapeutic classes of compounds. The calculated mean 
vector (i.e., the average vector representing a group of compounds with similar 

0 biological activity) is shown in Table 1 for hallucinogenic compounds. 
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TABLE 1 

Membrane binding properties or hallucinogens 
constitute a training set for analysis 



Compound 


1AM. PC 


IAM.PE 


1AM. PS 


IAM.SM 


n-acetvlmescaline 


0.706 


1.323 


0.736 


1.051 


* 

psilocin 


4.656 


4 818 


4.941 


7.532 


(R)-2,5-dimethoxy-4- 
iodoamphetamine mescaline 


55.53 


29.76 


260.200 


53.41 


2,5-dimethoxy-4- 
bromoamphetamine 


15.260 


20.950 


143.400 


20.970 


mescaline 


1.476 


2 134 


12.280 


2.074 


p-iodoamphetamine 


24.880 


33.660 


258.100 


32.410 


mean membrane binding 
properties (Mean vector) 


17.085 


5.441 

1 


113.276 


19.574 

i 



As shown in the last row of Table 1, the average vector representing 
the membrane binding properties of hallucinogenic compounds is {17.08 5 .441 

20 1 13.276 19.574} which corresponds to the membrane binding properties on 

{IAM PC, IAMPE, IAM.PS, IAM.SM}. The average vector of several other training 
(control) sets were calculated as shown for hallucinogens in Table 2. The goal is to 
compare the membrane binding constants of an unknown compound to the average 
vectors representing each class As shown in Table 2, when the angle between an 

25 unknown compound {i.e., 3,4-methylendioxyamphetamine (MDA) in this example) 

was calculated relative to the average vector for each therapeutic class of compounds, 
the closest match was hallucinogens and the closest receptor match was serotonin. 
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Comparison of MDA to the mean membrane binding 
properties of compounds within different therapeutic classes 



Activity 


Angle 


hallucinogens 


2 56 


antihypertensive 


2.770 


anorexic 


7.29 


stimulants 


15 27 


serotonin* 


23.26 


antidepressant 


27 77 


antiparkinson 


44 20 


anxiolytics 


48 16 


dopamine* 


49 33 


antipsychotic 


49 90 


y-amino butyric acid* 


57 87 


opioid 


60 46 


sedative-hypnotics 


61 89 


anticonvulsant 


63 13 



20 * receptor class. 

These matches indicate that the test compound may have biological properties similar 
to hallucinogenic compounds and may act at serotonin receptors, which is indeed the 
case of MDA 

25 In one embodiment of this invention membrane binding constants are measured 

by injecting the drugs into high performance liquid chromatographs (HPLC) using 
immobilized artificial membrane (IAM) columns using the following 1AM stationary 
phases: 

30 ete lAM PC CI0C3 , ^AM.PE 0100 , ^AM.PS 0100 , and ctlCT IAM.SM cl0C3 

The peak position of the drug in the chromatogram is a measure of the affinitv of the 

. ...... . . v.'inpounj ku i AM ^. 
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Capacity factors are proportional to equilibrium membrane binding constants K 
according to 



where 0 is the phase ratio of the column. Effectively, the set of average k' values for a 
group of compounds with similar biological function constitute a membrane affinity 
fingerprint (MAF M ) of the biological function. Thus MAF M denotes a mean vector 
whereby the vector components are average membrane binding constants for a group 
of compounds (a training set) with similar biological activity or functionality, e.g., oral 
absorption. 

Establishing Training Sets 

The term "training set," as used to describe and claim this invention, is a set of 
compounds of common biological properties used to define/calculate an MAF^ value 
for that biological property. Normally distributed membrane/membrane mimetic 
binding affinities in N-dimensional space occur when bivariate log plots of ak' (i.e. PC 
vs PE, PC vs PS, etc ) have an ellipsoidal shape. The ellipse shape indicates that the 
affinity data for a given set of compounds having a known clinical efficacy or function 
are normally distributed in 2 dimensional space In the case when a compound falls 
outside of the 0 95 quartile (3 standard deviations from the mean), it is considered to 
be an outlier Deletion of outliers increases the ellipticity of the bivariate data, thus 
improving the multivariate normality For this reason, the general procedure used to 
establish a training set of compounds is to continuously (1) remove outliers, and (2) 
redraw the ellipse plot until all compounds fall within the 0.95 quartile. Ellipse plots 
were thus used to remove outliers and confirm that the membrane binding data were 
normally distributed in N-dimensional space. Ellipse plots for all possible combinations 
of membrane binding constants were used to define each training set: PC vs PE, PC vs 



activity or function or clinical efficacy Typically the compounds in the training set 



k=0K 



PS, PC vs SM, PC vs PS, PE vs SM, and PS vs SM 




i compounds wmcn nave a mmmon t^oiogica; 
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composition are in a predetermined molar ratio. In one embodiment the control 
compounds in the training set composition are in a substantially equimolar ratio, 
optionally in combination with one or more external standards. Preferably the 
compounds in this training set selected so that the ellipse plot of the membrane mimetic 
5 binding data for the compounds in the training set composition are such that all 
compounds fall within the 0 95 quartile. The training set typically includes at least 
two, preferably at least three, most preferably 5 or more control compounds. The 
training set composition can be combined with test compounds (and optionally training 
set compositions for other biological properties) to prepare mixture for LC/MS 
10 analysis of membrane/membrane mimetic binding properties. The training set 

compositions in accordance with this invention provide ideal internal controls for direct 
comparison of test compound affinities with those exhibited by the compound 
members of the training set in LC/MS analytical protocols. 

15 Determining Mean Vectors 

The mean membrane binding affinities for each training set are calculated by 
averaging the PC, PE, PS, and SM as illustrated in Table 3 for selected hallucinogen 
compounds Multivariate Analysis of Variance (MANOVA) confirmed that greater 
than 90% of the mean membrane binding affinities of all therapeutic classes and 
20 mechanism of actions were different at the 0 05 confidence level. In other words, the 
mean vector for each training set developed for use in accordance with this invention 
was unique compared to all other training sets. 



Compound 


IAM.PC 


IAM.PE 


IAM.PS 


IAM.SM 


mescaline 


-0.477 


0.670 


2 43 


-0.715 


psilocin 


0 900 


1 79 


2.80 


0.502 


MDA 


0.298 


1.32 


3.28 


-0.0647 


DOB 


2.07 


3.22 


4.68 


1 62 



Table 3 Mean Vector or Membrane Affinity Fingerprint (MAF M ) for Hallucinogens 
The data in the above table is Ln weighted 
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Predicting the Activity of Unknown Compounds 

P 2 and 0 Values 

In addition to the vector analysis method described above for calculation of 
angle 0 values, D 2 and Q values can also be used for classifying/comparing compounds 
by vector analysis of binding data (and other quantitative physicochemical properties) 
Membrane binding data for an unknown compound is an observation vector y, and the 
mean membrane binding data for the i* training set of compounds is denoted as the 
vector y r The classification of an unknown compound involves calculating the mean 
squared distance, D 2 , between y and y r 

D 2 : - (y-yy S - 1 (y-yj 

where S, is the covariance matrix of the i* group of compounds. The unknown 
compound is assigned to the therapeutic class or receptor type with the smallest D 2 , 
value. 

An alternate method for classifying compounds is by use of Q values. 
Q,(y) = -V 2 In | SJ-'/i (y- y y S; x (y-y,, 

where | S,| is the determinant of the covariance matrix of the i* group. When Q values 
are used to classify compounds, the unknown compound is assigned to the therapeutic 
class with the highest Q value. 

Both D 2 and Q were used to classify the activity of each compound comprising 
all of the training sets A compound was considered to classify correctly when its 
known therapeutic class or mechanism of action was within the top two hits based on 
D 2 and Q Table 4 shows the D 2 and Q values for the classification of MDA (a 
hallucinogen) 
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Activity 






naitucinogen 




5 ~>1 


antiparkinson 


1 77 

1 / / 


-7 74 


anorexic 


Q 1 4 


""Xi . J o 


antihypertensive 


4 OQ 


-1 Q1 


anajgesic 


i n a 

1 \J o 




diuretic 


n q 

U 7 


7 


antineoplastic 


n qq i 




antidepressant 




1 0 7 


antipsychotic 


25.5 


-115 


nucleoside 


6 38 


-12.3 


antiviral 


31 7 


-13 5 


anticonvulsant 


52.0 


-24.0 


sedative hypnotic 


104 


-49 3 



Table 4. D 2 and Q values for MDA compared to all therapeutic classes. 

20 Order ofElution 

The relative order of elution is a measure of the relative affinity of each 
compound for each membrane mimetic surface. Unfortunately, theta, D 2 , and Q values 
for a compound of unknown activity do not necessarily reflect elution order. Thus a 
low theta and D 2 or large Q may occur even though a compound does not have the 

25 same relative order of elution exemplified by the mean vector of the most comparable 
training set. Consequently, in vivo activity can be predicted using not only theta, D 2 , 
and Q values, but also the LC elution order of the unknown compound on the 
respective solid phase substrates vs that of the mean vector for a known training set 



: a, dT,cn. . . . .t i u «. > ^ itcii pi euicting trie in \ i \o activ::v or a 
compound is the peak width associated with each membrane mimetic surface Peaks 
widths evaluate the on-off kinetics of drug membrane interactions For example the 
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unknown compound BDFA was evaluated for hallucinogenic activity. Its theta value 
relative to the hallucinogen MAJP was 3 .6 and its order of elution matched that of the 
hallucinogen MAF M However, the peak widths of BDFA on SM and PS are much 
larger than the corresponding peak width of MAF ^ In vivo studies have shown that 
5 BDFA is inactive Both the theta value and order of elution indicate that BDFA would 
be active in vivo Only the difference in the peak widths provided evidence that BDFA 
would be inactive in vivo. 

In accordance with another embodiment of the present invention there is 
provided a system for screening test compounds for probable biological properties. 

10 The system comprises two or more membrane mimetic surfaces, each having a unique 
amphiphilic composition. The system includes means for quantifying the interaction of 
test compounds and control compounds with each of the surfaces and assigning a 
numerical value characteristic of said quantified interaction of the test compound and 
each respective membrane surface. Typically, the screening system also includes a 

15 database comprising numerical values characteristic of the quantified interaction of 

selected control compounds or training sets of control compounds with the membrane 
surfaces. At least a portion of the selected control compounds have a predefined 
biological property The system includes an analyzer (e.g. a programmable computer) 
for comparing the numerical values determined for the test compounds with the 

20 corresponding numerical values for control compounds and identifying control 

compounds having numerical values best matching those for the test compound. In 
one embodiment the system includes a graphics algorithm for displaying the numerical 
values for the test compound and the numerical values for at least a portion of the 
control compounds as visibly distinguishable vector quantities. In another embodiment 

25 the analyzer includes an algorithm using vector calculus manipulation to identify the 
control compound or compounds having numerical value best matching those of the 
test compound Thus, for example, the system can be programmed to report all 
control compounds wherein © (see above) is less than 15°, more preferably less than 
1 0° . or where T> 2 U minim izer* n : ^ 1 ■ * ' 

d^tivitu's pvmhi! reiati\e mcm&iane affinities such that their membrane binding data 
form a readily discernible grouping when presented as vectors in "membrane space" 
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See Figs 1 and 2. Fig 1 shows vector plots in 3-dimensionaI "membrane space" for 
compounds which act at serotonin receptors: buspirone, cinanserin, mescaline, 
methysergide, plamphetamine, psilocin, quipazine, meterogline and others. Fig. 2 
shows vector plots in 3-dimensional "membrane space" for compounds that act at 
5 dopamine receptors: apomorphine, clozapine, domperidone, haloperidol, pridinol, 

SCH23390, SKF38393, spiperone, sulpiride, chlorpromazine, dihyrexidine, dopamine, 
perphenazine, prochlorperazine 

That phenomenon (the grouping of membrane mimetic binding data as 
vectors in "membrane space") forms the basis of the present invention and allows 

10 multiple membrane binding constants to serve as a predictor of biological activity 
While Figs. 1 and 2 are iiiustrated as the products of pattern matching by vector 
analysis, other art-recognized mathematical analytical techniques such as multivariant 
analysis and principal component analysis can be applied to membrane binding data to 
compare test compounds with control compounds with known biological function 

15 Vector analysis is particularly useful for pattern matching in accordance with this 

invention in that it can be carried out, albeit not graphically represented, in more than 
three dimensions. Such techniques can also be applied to data sets that include, in 
addition to data characteristic of multiple membrane affinities, other biologically 
significant molecular descriptors. Further, although membrane binding data acquisition 

20 is illustrated using immobilized artificial membrane chromatography supports in a high 
performance liquid chromatography system, other techniques can be used, for example, 
computer chips or similar devices with immobilized lipids, capillary zone 
electrophoresis columns coated with membrane lipids, Langmuir Blodget films, 
liposomes, and adsorbed monolayers of lipids on any surface, for example an AFM tip 

25 and evaluating the change in oscillation of the tip in the presence of test and control 
compounds. In addition to such empirical data acquisition techniques, numerical 
values characteristic of both MAF and non-MAF parameters can be obtained by 
computer calculations Thus data sets for use in accordance with this invention can 



the pjc>cni indention are numerical quantities that can be derived from a known 
chemical structure Examples include surface area, polar surface area, number of H- 
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bonds, topological indices, solubility, etc. These parameters alone do not predict the 
general membrane binding properties, i.e., MAFs, of compounds. The term 
"calculated MAF parameters" refers to a combination of calculated non-MAF 
parameters, with or without individual membrane binding constants, that can be used 
5 to predict Membrane Affinity Fingerprints (MAF) Since HPLC retention times can be 
used to calculate membrane binding constants, methods used to calculate retention 
times of solutes are actually calculating MAF parameters. An example of a calculated 
retention time for 12 compounds can be found in (Amie, D. Davidovic-Amie, D 
Trinajstic, N., J. Chem. Inf. Comput. 1995, 35, 136-139.) The term "simulated MAF" 

10 refers to a calculated MAF parameter obtained from Molecular Dynamics Simulations 
of compounds with IAMs, bilayer membranes, other membrane mimetic surfaces, or a 
force field characteristic of these membranes. 

In another embodiment of this invention there is provided novel 
carboxyl-functional, head group-protected phospholipids of the formula HOOC-W- 

15 OPO : OZ wherein Z is protected glyceryl, 2-(protected amino ethyl), 2-protected 
carboxyl-2-aminoethyl, 2-protected carboxyl-2-protected amino ethyl or a readily 
cleavable protecting group, and W comprises the lipid residue of a biologically 
significant compounds such as sterols, steroids, and fatty acids including, for example, 
retinoic acid, or W comprises the lipid residue of other amphiphilic molecules found in 

20 natural biological membranes including, but not limited to lecithins, lysolecithins, 

cephalins, sphingomyelin, cardiolipin, glycolipids, gangliosides, and cerebrosides Such 
compounds are useful for preparation of certain of the immobilized artificial 
membranes useful for comparing membrane interactions and predicting biological 
activities of test compounds Alternatively W is the lipid residue of a phospholipid of 

25 the general formula W-OP0 2 OB, wherein W can be acylglyceryl, diacylglyceryl, or N- 
acyl 3-0-(protected) sphingosin-l-yl, wherein "acyl" is C t -C 24 alkanoyl or C 8 -C 24 
alkenoyl 

In another aspect of this invention novel carboxyl-functional, head 

HOOC -V\ -OPO ; OCHX H : N (CH v h in the presence of an excess of a protected 
alcohol of the formula ZOH wherein Z is protected glyceryl, 2-(protected amino)ethyl. 
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2-protected carboxyl-2-amino ethyl or the residue of an acid protecting group, for 
example, 4-nitrobenzyl or 4-nitrophenylethyl. W is a lipid reside as defined above with 
the proviso that it is selected so that the starting compounds serve as a substrate for 
phospholipase D activity. Typically about a 3-fold to a 15-fold stoichiometric excess 
5 of the alcohol is used, and the reaction is carried out in a buffered aqueous organic 
solvent, for example chloroform or ethyl acetate, in the presence of a water soluble 
calcium salt The product carboxyl-functional compounds are useful for preparing 
novel head group-protected immobilized artificial membranes on surfaces of solid 
substrates having carboxyl reactive functional groups such as hydroxy, amino, and thiol 

10 groups. Removal of the protecting groups on the immobilized phospholipids provide 
the corresponding deprotected artificial membrane structures. 

The term ''protected" or "protecting group" as used in describing the 
present invention refers to those chemical moieties that 1) can be temporarily bonded 
to a reactive functional group {e.g., phospho, carboxyl, hydroxy or amino) to prevent 

15 subsequent unwanted reaction of the reactive functional group during other chemical 
modification of the compound bearing said functional group and that 2) can be 
removed from the functional group or groups at a subsequent synthesis step without 
unwanted reaction of other portions of the protected molecule Protecting groups for 
amino, carboxyl and hydroxy functionalities are well known in the art, as are their 

20 respective synthesis and conditions for removal (deprotection). Preferred protected 
alcohols useful for the preparation of the carboxyl functional, head group-protected 
phospholipids in accordance with this invention are ispropylidene glycerol, allyl serine, 
2-(t-butoxycarbonylamino)ethyl, 4-nitrobenzyl alcohol, and 2-(4-nitrophenyl)ethanol 
Another artificial membrane stationary phase useful for assessing 

25 membrane interaction and predicting biological properties in accordance with this 

invention is a novel immobilized ceramide-based membrane In one embodiment the 
ceramide stationary phase is prepared by covalently bonding N-(13- 
carboxyltridecanoylVD- erythro -sphingosine through the ovcarboxyl group (via the 
ceramide-imidazolide) to silica propvlamine and suhseauent Ci ^nrl r ^ o^-^n.-/ 
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The immobilized artificial membranes used in accordance with this 
invention to predict biological activities can be prepared using one or more than one 
o>-carboxyl substituted lipids. Artificial membrane structures useful as a stationary 
phase for HPLC chromatographic determination of membrane interaction and having a 
predefined ratio or stoichiometry of phospholipid components can be prepared to 
simulate those lipid/phospholipid ratios known to exist in a predetermined biological 
membrane Thus, the binding affinities of test compounds and control compounds to 
both homologous and heterologous membrane substrates can be measured and 
compared to assess/predict biological function. 

Membrane binding data (numerical values relating to interaction with 
membranes or membrane mimetic surfaces) were obtained using HPLC 
chromatography on a multiplicity of membrane mimetic surfaces, preferably 
immobilized artificial membranes (IAMs). The membrane mimetic surfaces used to 
quantitate membrane lipid interactions were prepared from analogs of 
phosphatidylcholine (PC), phosphatidylserine (PS), phosphatidylethanolamine (PE), 
and sphingomyelin (SM). These membrane lipids were chosen because of their 
established asymmetry in plasma membranes. Synthetic methods for preparing 
immobilized membranes are well established. Experimentally, drugs were injected into 
high performance liquid chromatography (HPLC) columns. The peak positions in the 
chromatograms measure the affinity of compounds for the immobilized membrane lipid 
surfaces. 

Comparing membrane binding properties among a multiplicity of 
compounds required a pattern matching process whereby the patterns being matched 
were sets of membrane binding data. Compounds with similar pharmacology have 
been found to cluster in membrane space. This supports the hypothesis that 
compounds with similar pharmacology have similar membrane binding properties. The 
dispersion of vectors among the compounds of a therapeutic class defines the region of 
membrane space where activity of that therapeutic class resides When outliers in 
membrane space are observed for member of the class, other factors in addition to 

. ,. , llf v ^ ,o . - . ^MM>ier,i v\ \\ n t ne>c ^ omniums ^ exniDiting ef.senr.u. 

membrane binding properties needed for drug penetration through the blood brain 
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barrier Although active transport and efflux mechanisms have recently been 
emphasized as regulating penetration through the CNS, membrane binding must 
significantly contribute to the regulation of CNS transport because all CNS compounds 
tested have unique membrane binding compared to non-CNS compounds. The region 
5 occupied by non-CNS compounds differs from that of CNS compounds both in vector 
direction and vector magnitude. Thus, sets of membrane binding data can be used to 
determine if a given compound is likely to penetrate the blood brain barrier. Once 
CNS compounds penetrate the blood brain barrier and are free to distribute into 
different regions of the brain, membrane binding properties unique to each therapeutic 

10 class regulates the localization of compounds within the CNS. Specifically, when 
compound localization within the CNS is regulated primarily by membrane binding 
properties, then drug vectors within therapeutic classes are expected to occupy a 
unique region of membrane space. This clustering was observed for all therapeutic 
classes of compounds tested which suggests that the in vitro membrane binding model 

15 has significant utility in modeling the distribution of compounds within the CNS. 

Vector plots in membrane space are useful for demonstrating that the 
angle between two drug vectors is small when compounds exhibit similar biological 
activities However, angle calculations are not restricted to three membrane binding 
constants. Better grouping by therapeutic activity is attainable when more than three 

20 membrane binding constants are pattern matched. For instance, diazepam pattern 

matches membrane binding properties of anxiolytic compounds much better when four 
membrane binding constants are used for the analysis. As shown in Table 5, when only 
three membrane binding parameters were used for pattern matching diazepam to the 
database, only two anxiolytics were among the top ten compounds, whereas five 

25 anxiolytics were identified when four membrane binding parameters were used 

Although larger angles are calculated when four verses three membrane parameters are 
used (Table 5), pattern matching membrane binding constants in N-dimensional 
membrane space appears to give better groupings of compounds regarding biological 
activity The observation that compounds with bioloeical activities HifTprpr-t 

, v\cci. jiiiereni ^ V> tnerapeutic wx i> unavoidable in addition, trie iack oi a 
non CNS compound among the top ten matches in Table 5 emphasizes the clear 
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distinction that exists between the membrane space occupied by CNS and non-CNS 
compounds. 



5 Table 5 

Pattern Matching Membrane Binding Constants in 3 Dimensions and 4 Dimensions. 
Angles were calculated from membrane binding data and are relative to the membrane 
binding data of diazepam. The Membrane Space for each analysis is given. 



10 



Three Dimensional Pattern Matching 


Four Dimensional Pattern Matching 


Membrane Space = 


{PC, PE, PS} 


Membrane Space = 


{PC, PE, PS, SM} 


compound 


angle 


activity 


activity 


angle 


compound 


diazepam 


0 


anxiolytic 


anxiolytic 


0 


diazepam 


(ref. compound) 










(ref. compound) 


promethazine 


1 48 


depressant 


anticonvulsants 


2.83 


nitrazepam 


chl orodi azepoxide 


1.82 


anxiolytic 


anxiolytic 


3.34 


chl orodi azepoxide 


2Br-LSD 


2.56 


LSD antagonist 


LSD antagonist 


5.71 


2Br-LSD 


nitrazepam 


2.45 


anticonvulsant 


anxiolytic 


5.77 


oxepam 


procyclidine 


3.04 


nonCNS 


antidepressant 


5.83 


R-deprenyl 


p razepam 


3.16 


anxiolytic 


anxiolytic 


6.48 


medzepam 


cyclobenzaprine 


3.28 


psychosedative 


hallucinogen 


7.38 


psilocin 


imipramine 


3.44 


antidepressant 


anxiolytic 


8.41 


prazepam 


nefopam 


3.78 


antidepressant 


anxiolytic 


10.15 


halazepam 
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Synthesis of Ligands for Preparing Immobilized Artificial Membranes 

l-myristoyl-2-(l 3 -carboxytridecanoyI)-5«-3 -phosphatidylcholine (PC- 
30 COOH) was used as a common chemical intermediate to prepare immobilized artificial 
membrane (IAM) surfaces for measuring drug-membrane interactions. Previous 
preparation of IAMs containing PC, PG, PA, PE, and PS from ether phospholipids 
required - 30 steps, whereas, the use of PC-COOH to enzymatically prepare protected 
o>-carboxyl phospholipid ligands (pl 2 ^ 00 ^ from phospholipase D (PLD) reduced 
35 the total number of synthetic steps to ~ 20. The five pl u>-cooh1> ligands prepared were 
analogs of phosphatidylglycerol (PG 2fa > COOHP , la), phosphatidylserine (ps 2 ^ 00 " 7 lb). 



40 



^ . a. , ii^pnaiiu\ lation reactions usinn H*. - 

COOH were typically - 90-95° o complete An interesting finding was that serine allyl 
ester was an excellent substrate for PLD and quantitative transphosphatidylation of 
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PC-COOH with serine allyl ester was obtained. This is the first report describing the 
preparation of SM analogs suitable for bonding to molecular surfaces. PL"" 000 *^ 
protecting groups were removed after bonding the ligands to silica propylamine and 
endcapping with straight chain anhydrides After deprotection, five immobilized 
5 membrane surfaces were prepared that contain interfacial functional groups identical to 
PG, PA, PE, PS, and SM membrane lipids. 

Although single chain ether PLs were previously used to prepare the 
first generation of I AM surfaces, two factors prohibited the wide spread use of ether 
lipids to prepare IAMs for measuring drug-membrane binding constants. One problem 

10 was that optimum predictions required immobilized lipids that are structurally similar 
to endogenous lipids found in cell membranes and single chain ether lipids lack several 
interfacial functional groups commonly found in the majority of membrane 
phospholipids. The second problem was that ether PLs required numerous synthetic 
steps and the preparation costs are high. Accordingly, in another embodiment of the 

15 invention efficient synthetic strategies were developed to prepare carboxyl analogs of 
membrane phospholipids for immobilization that bear identical interfacial functional 
groups as those found in the PLs forming the cell membranes. 
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interfacial polar groups 



, t 




^ ^COOAII 

^^-"^ P«i 2tD-COOH/P 

(1d) I ^ b 

NHBoc 



Scheme l Structures of protected diacylated phospholipid analogs ofPG, PA, PE 
and PS After immobilization and following deprotection, these ester PLs have 
identical interfacial functional groups to endogenous membrane phospholipids (PG, 
10 PA, PE, or PS) 
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Scheme 1 shows the protected diacylated pl 2uM::oohp ligands prepared 
during this work The protecting groups, Rp, shown in Scheme 1 are the same 
protecting groups used to prepare single chain ether PLs that were immobilized except 
for the serine analog Ether PLs containing serine carboxyl were protected with a tert- 
5 butyl ester, whereas the present synthetic procedure used an allyl ester group because 
this protecting group does not interfere with PLD transphosphatidylation reactions. 
Scheme 2 shows the carboxyl SM analog (2) prepared as well as the main structural 
differences between SM^°° HP and vV^°° Hr? ligands The SM^°° HT ligand was 
used to prepare an IAM surface that model drug binding to endogenous membranes 
1 0 enriched in SM, an example is the kidney brush border cells where 80% of SM in the 
plasma membrane is located in the outside membrane leaflet. 
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The general synthetic route to prepare double chain ester PLs, instead 
of single chain ether PLs, bearing u>-carboxyl groups is shown in Scheme 3. PLA 2 
hydrolysis of DMPC to obtain lysolecithin (5) was found to be quantitative but 
reaction conditions are enzyme specific In addition, PLA 2 s from different 
5 manufacturing lots contain different contaminants that interfere with product 

formation. Several PLA 2 s including Naja mocambique mocambique PLA 2 (pi = 8.8), 
Naja mocambique mocambique PLA 2 acidic isoenzyme and bee venom PLA 2 (Apis 
mellifera) were tested for this work and all enzymes could be used to obtain 
quantitative hydrolysis of DMPC, but, reaction conditions needed to be optimized for 

10 each enzyme Bee venom PLA 2 was reported to have high activity in low water- 
organic solvent systems (1 7% buffer/chloroform V/V). Preliminary studies were 
intended to determine if 1 .7 % V/V buffer/CHCl 3 was sufficient for scale up reactions 
whereby several gram quantities of 5 could be produced from a single reaction. The 
first reaction was a 100 fold scale up from the published procedure and Tris buffer was 

15 used instead of HEPES. DMPC (1.0 g) incubated with bee venom PLA 2 (1000 units) 
in 1.7% V/V Tris buffer (50mM, pH 7.2)/CHCl 3 did not react over 4 hours, and when 
the aqueous phase was adjusted to 5% V/V buffer/CHCl 3 the reaction still did not 
proceed Only when 10% V/V buffer/CHCl 3 conditions were used did the reaction 
proceed to completion in - 6 hours. Product precipitation in acetone always resulted 

20 in large amounts of the buffer reagent contaminating 5. For this reason, reactions were 
investigated to identify an acetone soluble buffer that does not interfere with PLA 2 
cleavage of phospholipids. Triethanolamine buffer was found to be acetone soluble 
and permitted PLA 2 hydrolysis of DMPC reactions to proceed quantitatively to 
generate lysolecithin (5). It was found that in triethanolamine buffer, Naja 

25 mocambique mocambique PLA 2 (pi = 8 8) had the highest efficiency for DMPC 

hydrolysis among the enzymes tested which included Naja mocambique mocambique 
PLA 2 acidic isoenzyme and bee venom PLA 2 (Apis mellifera) Lysolecithin purification 
required three acetone precipitations Several PCs were also hydrolyzed in 
quantitative vields usinc PI A- indicate r r ^ ^ >^ ' j - x 1 

pr-CQQH (4, 1 v\as prepared horn 5 as described Larher studies 
indicated that PLD from either vegetables or microbacterials has activity in 
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organic/aqueous systems and it accepts unnatural substrates (including both primary 
and secondary alcohols) and generates phosphoester products. Therefore we extended 
PLD transphosphatidylation reactions to prepare the final ligands la-Id (Scheme 1). 
Ligands la and lb were prepared directly from PLD and did not require additional 
5 headgroup protection steps prior to ligand immobilization In contrast, ligands lc and 
Id require additional synthetic steps to protect the primary amines on each of these 
compounds so that these amines do not interfere with immobilization of the pl 2u> ^ oohp 
ligands. The major side product in PLD reactions is the formation of PA because 
water is an efficient competitor to the alcohols used for the PLD reactions. Thus one 
10 of the primary goals during PLD reactions is to minimize formation of PA side 

products regardless of PC anaiog subject to PLD transphosphatidylation. Based on 
TLC, PLD reactions using compound 4 show either no or trace amounts of PA- 
COOH 

Commercially available IPG alcohol was used to prepare PG 2a ^°° HT> 
15 (la) by PLD transphosphatidylation; the reaction was facile and an -83.2% yield was 
obtained after product purification by flash chromatography. The PLD reaction was > 
90 % complete and product recovery was decreased due to column chromatography. 
No PA side product formed during the reaction conditions used for PLD conversion of 
compound 4 to the ligand la. Although compound 4 is not very soluble in buffer/ethyl 
20 acetate solvents, the yield of PG 2w * COOH/p was much higher and the reaction times was 
much shorter, in this mixed solvent system compared to buffer/CHCl 3 solvent system 
Similar observations were found for the PLD transphosphatidylation reactions using 
compound 4 and ethanolamine which was the alcohol used to prepare pe 2 ^ 00 ^ (i c ) 
The most important consideration when using PLD enzymology to 
25 prepare ps 2 "^ 00 * 17 (id) is the presence of both the serine carboxyl group and the u>- 
carboxyl group located on the terminal part of the acyl chain which is used for 
immobilization of the ligand The serine carboxyl group must be protected prior to 
PLD reactions so that the u>-carboxyl remains free and can be used for immobilization 

. ■ ' ' ■■ - - ^ ^ ■ ( _ _ ^ , , . , . ^ _ , i l : i ] i : ; n ; j l $ u r I a l e t i ■ j. 1 

. c-^nonj* using serine meth\l ester, cannot be used because conditions for deprotecting 
the serine methyl ester will also hydrolyze the interfacial esters of the immobilized 
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lipids 

Thus a specific protecting group for the serine carboxyl had to be identified 
that permits not only the PLD transphosphatidylation to proceed efficiently, but, the 
protecting group must be able to be removed under mild deprotection conditions to 
5 generate the free serine caxboxyl after immobilization. Thus synthetic strategies for 
preparing the ligand Id by PLD enzymology tested serine benzyl ester, phenacyl ester 
and allyl ester All of these esters can be removed under mild solution conditions 
Steric hindrance of benzyl alcohol and 2-hydroxyacetophenone caused these alcohols 
to be poor substrates for PLD. However, serine allyl ester underwent quantitative 

10 PLD transphosphatidylation and this chemical intermediate was not purified, the 

unreacted serine amine was directly proiecied with Boc to form the ligand Id which is 
suitable for immobilization After immobilization, allyl ester can be removed using 
tributyltin hydride (Bu 3 SnH) and the catalyst bis(triphenylphosphine)palladium (II) 
chloride PdCl 2 (PPH 3 ) 3 The one pot two step reaction from compound 4 to ps 2tl>COOHP 

15 (Id) had a overall yield of 87.1% after the product was purified by column 

chromatography. The reason for not protecting the serine amine prior to PLD reaction 
is that a free amino is usually needed for serine to be a good substrate for PLD 
enzymes 

The preparation of SM^ 00 ^ (2) (see Scheme 4) ligand involved 
20 aminolysis of endogenous SM to form sphingosine-l-phosphocholine with a free 

amine Prior to reacting the free amine with a large ring cyclic anhydride, the allylic 
alcohol was protected with ter/-butyldiphenylsilane (TBDPS) The TBDPS protecting 
groups can be removed with TBAF which does not degrade the silica surface after the 
ligand 2 is bonded to silica 
25 The ester ligands were bonded to silica surfaces similar to procedures 

previously developed in our laboratory See U.S. Patent No. 4,931,498. Typically 
IAMs are prepared from 3 sequential bonding processes. The first bonding process 
links the PL*^ 00 ^ ligands to the surface using CDI as the activation agent and FTIR 
^pectro^cnpv U " ^p- 1 ♦ - ' r - * ■ ■ 
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SCHEME 4 
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A detailed study of endcapping reactions has been performed, and the 
results from this and other work indicate that either long chain or short chain 
anhydrides can be used for endcapping However, a critical concept is that the 
endcapping reactions are performed until the surfaces are ninhydrin negative 
5 Typically, decanoic (10) anhydride is the first endcapping reagent followed by 

propionic (C3) anhydride. However, double endcapping with C3 or CIO or any other 
combination can be performed depending on the use of the IAM packing materials. 
For instance, NMR studies of the IAMs require that the IAMs suspend in D 2 0 which 
can only be achieved by omitting the endcapping with CIO anhydrides. Thus 

10 endcapping is not routine from a synthetic perspective because it can be used to 
control the hydrophobicity, wettability, and other properties of the final surfaces 

Since both endcapping may be varied and the protection groups vary 
among pl 2 ^°° hp ligands, a nomenclature is needed to identify the chemical surface 
intermediates formed throughout the multistep bonding process. PL 2 *^ 00 " 71 * ligands 

15 are always bonded to the silica surface first, followed by endcapping and deprotection. 
The nomenclature throughout the sequential bonding processes for preparing the 
^AM.PS 03 C3 surface (10) was as follows: 

SPA -> c,te TAM.PS p -> ttlCT IAM.PS p/C3 -> CTlCT IAM.PS p/C3C3 -> ^AMPS 0303 (10) 

20 

where SPA is silica propylamine starting material which becomes •* tCT IAM PS p after 
immobilizing ps 2 "" 000 "* which then becomes ^AM.PS* 03 ' 03 after two sequential 
endcapping steps using C3 anhydride The superscript p indicates that serine 
protecting groups have not been removed After removal of all the protecting groups, 
25 e * lCT IAM PS P C3C3 becomes the final product 10. This nomenclature, which is 

exemplified for the preparation of e * la lAM.PS C3 C3 surfaces, is used for all IAMs 

Bonding strategies of the pl 2 ^ 00 * 1 * ligands to the silica surfaces were 
the same as that of PC-COOH which did not require protecting groups However, the 
pl **coohp i igands are ] eS s soluble in CHC1 than pr-fOOH - 

,,iitd prop\idmme was lacne U hen ether PLs were used for preparing IAMs 
aqueous HCI was used to deprotect Boc, tert-but\\ esters, and isopropylidene (IPG) 
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groups because the immobilized ether lipids were stable to this reaction condition. 
However, the ester PLs described in this report are susceptible to acid hydrolysis and 
we used anhydrous TFA to remove acid sensitive protecting groups. 

We have found that the only way to monitor trace amounts of 
5 functional groups on the I AM surface after the pl 2 ^ 00 *^ ligands are immobilized is 
to use ER microscopy. However compared to all other protecting groups tested, direct 
monitoring the deprotection of the allyl ester groups on e * ICT IAMPS PC3 /C3 by IR 
microscopy was impossible because allyl ester group does not have an IR band in a 
unique region of IR spectrum. From peptide chemistry, deprotection of allylic groups 
10 are quantitative using tributyltin hydride. Nevertheless the deprotection was essential 
to monitor because interfacial surface chemistry significantly affects many chemicai 
reactions and deprotection of allyl ester may not be facile on IAM surfaces. 
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IR band shifts with serine carboxyl group on ^AM.PS 0303 in the ER 
spectra under acidic and basic conditions were used to monitor the extent of 
deprotection Thus after deprotecting ^AM PS P C3 C3 with Bu 3 SnH and PdCl 2 (PPh 3 ) 3 , 
IR spectra of the acidified I AM surface (serine-COOH) and basified I AM surface 
5 (serine-COO") showed a decrease in the integrated ester region between 1780-1680 
cm' 1 by - 1 .3 fold. This decrease in intensity was expected because the surface 
contains two interfacial esters from the immobilized phospholipid. The serine carboxyl 
shift from 1710 cm' 1 to 1550 cm 1 when the serine COOH is ionized to serine COO* 
chemical species moved the deprotected functional group to a new region of the IR 

10 spectra allowing the deprotection reaction to be monitored. 

Numerous synthetic strategies can be developed to prepare the 
assortment of ester pl 2uhCOOHP ligands shown in Scheme 1 For instance, PLA 2 
hydrolysis of endogenous phospholipids including PG, PA, PE, PS can be performed 
followed by protecting polar headgroups and linking a long chain acyl group with an 

1 5 co-carboxyl group. This type of synthetic strategy whereby each p\J^°° nv Hgand is 
synthesized from an endogenous precursor has several disadvantages, the most 
obvious of which is the generation of several extra synthetic steps to obtain all PL 2t ^ 
cooh? Hgands 

However, equally important considerations are (1) endogenous starting 
20 materials like PG, PA, PE, and PS are expensive, and (2) the lysophospholipids of PG, 
PA, PE, and PS have poor solubility in aprotic solvents which makes them difficult to 
undergo following synthetic chemical reactions. The selection of PC-COOH (4) as a 
common intermediate (Scheme 3) to prepare ester pl 2 ^ 00 ^ thus avoids the usage of 
expensive starting materials PG, PA, PE, and PS and eliminates lysophoshpolipid 
25 intermediates The synthetic strategy described herein is efficient for both small and 
large scale reactions and will be generally useful for preparing membrane surfaces. 
The ligands can be bonded to any molecular surface 

General Experimental Procedures 

^'oufczht trnm i Tcnzyrne Pharmaceuticals ; L arribndgc, MA; Bee venom PLA : ^Api> 
mellifera), Naja mocambique mocambique PLA : (pi = 8 8), Naja mocambique 
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mocambique PLA 2 acidic isozyme, PLD (from Streptomyces species), CaCk 2H 2 0 and 
ethanolamine were ordered from Sigma Chemical Co. (St. Louis, MO) 
Triethanolamine (TEA), tris(hydroxymethyl)aminomethane (Tris), 
ethylenediaminetetraacetic acid disodium salt (EDTA), anhydrous trifluoroacetic acid 
5 (TFA), 1,-1-carbonyldiimidazole (CDI), propionic anhydride (C3), decanoic anhydride 
(CIO), tributyltin hydride (Bu 3 SnH), bis(triphenylphosphine)palladium (II) chloride 
(PdCl 2 (PPh 3 ) 3 ), alcohol-free anhydrous chloroform (stabilized with amylene), 
methylene chloride, 4-dimethylaminopyridine (DMAP), sodium acetate, di-ter/-butyl 
dicarbonate (Boc anhydride), p-nitrophenethyl alcohol (NPEA), triethylamine, tert- 

10 butyldiphenylsilyl chloride (TBDPSC1), imidazole, w-butanol, 1, 12- 

dodecanedicarboxylic acid, tetrabutylammonium fluoride (TBAF) and aiiyi bromide 
were purchased from Aldrich (Milwaukee, WI). Sphingomyelin (brain) and eggPC 
were ordered from Avanti Polar Lipids, Inc. (Alabaster, AL). Deuterated solvents for 
NMR spectroscopy including methanol and chloroform were from Cambridge Isotope 

15 Laboratories, Inc .(Andover, MA) Chloroform (alcohol free), acetone, 

isopropylalcohol, methanol, acetone, hexane, ethyl acetate, acetic acid, hydrochloride 
and sodium hydroxide were from Mallincrodt, Inc (Paris, NY). Labmotor, stirring 
shaft and stirrer blades were purchased from ACE Glass, Inc (Louisville, KY). 
Phospray reagent used for detecting phosphorus-containing organic compounds and 

20 ninhydrin reagent used for detecting amine-containing compounds were ordered from 
Supelco, Inc (Bellefonte, PA). Silica propylamine (SPA) was provided by Rockland 
Technologies Inc (Newport, DE) SPA (RN39-94) was used to prepare ""lAM.PL 
surfaces and it is trifunctional and has a surface area of 180m 2 /g and a pore size of 80 
Angstroms The density of propylamine groups on RN39-94 is 3 09 jimol/m 2 SPA 

25 (RN38-94) was used to prepare IAM SM surfaces and it is monofunctional and has the 

same surface area and pore size as RN39-94. The density of propylamine groups on 

RN38-94 is 2 06 (ig/nr All infrared spectrum were recorded on a Nicolet magna-ER 

spectrometer interfaced with a IR Plan I microscope Elemental analysis was 
performed on a Perkin-EImer PF 240 in the \firroamK-+^-^i t -aK - 
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l-myristoyl-s/i-2-hydroxyphosphatidylchoIine (5, lysolecithin) 

DMPC (10 0 g) was weighed and dissolved in chloroform (300 ml) 
Then 32.0 ml of 50 mM TEA (pH 7.2) buffer containing 2.5 mM Ca~ and 0.25 mM 
EDTA was added. Naja mocambique mocambique PLA 2 enzyme solution (1 ml, 1,000 
5 units/ml in 50 mM TEA buffer) was pipetted into the DMPC suspension. The viscous 
white reaction was paddle-stirred at 100 rpm for either 7.5 hours (basic isozyme) or 
1 1 .0 hours (acidic isozyme). Isopropylalcohol (100 ml) was added to the reaction 
mixture before rotoevaporation to avoid bubbling during solvent removal. After 
rotoevaporation, the residue was suspended in methanol (100 ml) and sonicated (1-2 

10 min) to solubilize the residue. The methanol solution containing product was slowly 
pipetted into acetone (- 500 ml) which caused a white precipitate (the product) to 
immediately form 

Complete precipitation required overnight incubation at 4 °C. The 
precipitated product was collected by filtration using an F glass-sintered filter. TLC 

15 (CHCl 3 /CH 3 OH/H 2 0 = 65:35:4) showed trace amounts of myristic acid and TEA were 
contaminating the lysolecithin product and a second precipitation was needed to obtain 
5 that did not contain any traces of myristic acid A third precipitation with the pH 
adjusted to 9 7 completely removed the TEA from the product and pure compound 5 
(5.16 g, 73 .1%) was obtained after vacuum dry overnight. Recovery of lysolecithin 

20 from the 3rd filtrate by an additional precipitation increased the yield to -80.0%. This 
reaction was performed ~ 50 times to obtain 285 g of lysolecithin which was used to 
prepare PC-COOH (4) by an established procedure. TLC in CHC1 3 /CH 3 0H/H 2 0 
(65 25 4) showed a pure product (Rf = 0 34); FTIR (thin film, cm' 1 ) 2955.24, 
2917 03,2849 68, 1734 22, 1467.66, 1234 61, 1086 39, 1056 77. l H NMR (300 

25 MHz, CDC1 3 , 6 ppm): FAB-MS [M+tT] calcd for C^H^PN 468.6, found 468 5. 

l-Myristoyl-2-(13-carboxytridecanoyI)-5«-3-phosphatidyIcholine (4, PC-COOH) 

L-serine ally! ester (3) 

■-leThviamine 1 2^ mi > were addec -\rur stn i ing ior 24 hours, the reaction mixture was 
filtered, and the filtrate was concentrated via rotoevaporation to yield a yellow oil 
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The crude product was purified by flash chromatography using a step hexane/ethyl 
acetate gradient from hexane (500 ml) followed by hexane/ethyl acetate (4:1, 1000 ml), 
hexane/ethyl acetate (2:1, 1000 ml), and hexane/ethyl acetate (1:1, 1000 ml). After 
drying in a vacuum oven overnight, 27.2 g of the intermediate N-ferf-butoxycarbonyl- 
5 L-serine allyl ester was obtained (93% yield). TLC in ethyl acetate/hexane (1:1) 

showed a pure product (Rf= 0.39); FTIR: (thin film, cm' 1 ) 3415.86, 2978.14, 2933.60, 
2886.66, 1743.43, 1715.75, 1512.33, 1392.96, 1370.85, 1342.40, 1165.13, 1060.54. 
l H NMR (300 MHz, CDC1 3 , 6 ppm): 5.87 (1H, m), 5.54 (1H, d, J - 7.79 Hz), 5.28 
(1H, d, J= 17 17 Hz), 5.20 (1H, d, J = 10 40 Hz), 4.61 (2H, d, J = 5.62 Hz), 3.87 

10 (2H, m), 1 40 (9H, s) CI-MS [M+H + ] calcd for C U H 19 N0 5 246.3, found 246.0. 

N-terf-butoxycarbonyl-L-serine ally! ester (27.2 g) was placed in a 500 
ml r. b. flask and cooled to 0 °C To this yellow oil, 50% TFA/CH 2 C1 2 (100 ml) was 
added under positive N 2 pressure. The reaction mixture was stirred at 0 °C (30 min), 
warmed to room temperature, and stirred for 2 h. TFA and CH 2 C1 2 were removed by 

15 rotoevaporation to yield a yellow oil. The yellow oil was dissolved in ether (200 ml). 
After 3 h at 4 °C, the product completely precipitated. The precipitated product was 
collected by filtration and washed with ether (50 ml) and hexane (20 ml). The white 
precipitate was dried in a vacuum oven overnight. Dry L-serine allyl ester (3) was 
obtained (15.6 g, 100% yield). TLC in CH 2 C1 2 /CH 3 0H/H 2 0 (65:25:4) showed a pure 

20 product (Rf = 0.58). FTIR: (thin film, cm" 1 ) 2994.54, 2960.64, 2891.49, 1747.28, 

1674 06, 1201.52, 11.45.92. J H NMR (300 MHz, CDC1 3 , 6 ppm): 5.86 (1H, m), 5.30 
(2H, m), 4 59 (2H, d, J = 6.43 Hz), 4.09 (2H, m), 3.77 (1H, m). FAB-MS [M+H + ] 
calcd for C 6 H n N0 3 146.2, found 146.0. 

25 3-0-(tert-but>ldiphenylsiIyI)-sphingosine-l-phosphocholine (6) 

Sphingosine-l-phosphocholine was prepared by hydrolysis of SM with 
minor modification Briefly, a solution of SM (10.5 g, 15 mmol), «-BuOH (150 ml), 
and HC1 (45 ml, 6 0 N) was stirred at 100 °C for 1 hour followed by rotoevaporation 
and vacuum drvine The sphincosine- 1 -phosphnrhniinp nrnrf* • ^ -i; - ♦ i- * 

nospnochoime wa^ rh^niv ed m I f T ? ( ^ [ iM mil followed bv tue audi li on of TBDPSCi 
(23 4 ml. 90 mmol) and imidazole (12 2 g, 90 mmol) in that order The reaction was 
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stirred (40 °C, 10 h) followed by rotoevaporation and purification by column flash 
chromatography. A white solid (5 64, overall yield 53%) was obtained after column 
chromatography TLC in CHC1 3 /CH 3 0H/H 2 0 (65:25:4) showed a pure product (Rf = 
0 25) FTIR: (thin film, cm 1 ) 3390.56, 3071.58, 3046 65, 3015.28, 2956.99, 2925.18, 
5 2854 .48, 1963.80, 1898.78, 1830.79, 1620 98, 1428.53, 1239.36, 1059.06, 973.34. 
*H NMR (300 MHz, CDC1„ 6 ppm): 7 66 (4H, m), 7.32 (6H, m), 5.43-5.09 (2H, m), 
4.45 (3H, m), 4 12 (4H, m), 3 40 (9H, s), 3.39 (1H, m), 1.64 (2H, m), 1.24 (22H, m, 
br s), 1 .00 (9H, s), 0.86 (3H, tr, J = 6.45 Hz). FAB-MS [M+H + ] calcd for 
Cj^NX^PSi 703 1, found 703.2 

10 

3-0-(tert-buryldiphenylsilyl)-sphingomyelin (2, SM"^ 00 "*) 

A 200 ml solution of 6 (5 g, 7.1 mmol) and 1,12-dodecanedicarboxylic 
anhydride (prepared from 1,12-dodecanedicarboxylic acid (5 0 g, 19.35 mmol) and 
DCC (4 0 g, 19 39 mmol) as described [Pidgeon, 1989 #8]) was stirred at r.t. for 10 

15 hours The reaction solvent was removed by rotoevaporation and product was purified 
using a step gradient mobile phases consisting of CH 2 C1 2 (500 ml), CH 2 C1 2 /THF (1:1, 
500 ml), CH 2 Cl 2 /CH 3 OH (1:1, 1000 ml), and CH 2 C1 2 7CH 3 0H/H 2 0 (65:35:4, 3000 ml). 
The chromatography fractions containing the product were pooled, and after 
rotoevaporation and vacuum oven drying 5 4 g of SM" 000 " 1, (2) was obtained (80.6 

20 % yield) TLC in CHC1 3 /CH 3 0H/H 2 0 (65:25:4) showed a pure product (Rf = 0.47) 
FTIR: (thin film, cm" 1 ) 33 12.05, 3071 70, 3046.10, 2923.61, 2855.04, 1960.21, 
1899.87, 1831.31, 1742.70, 1703.39, 1660.42, 1546.28, 1466.74, 1428.35, 1065.80, 
967.63 'H NMR (500 MHz, CDC1 3 , 6 ppm): 7 64 (4H, m), 7 34 (6H, m), 5.34-5 04 
(2H, m), 4.39 (1H, dd), 4 24 (2H, m), 4 02 (2H, m), 3 68 (3H, m), 3.21 (9H, s), 2 30- 

25 1.92 (4H, m), 1 76-1.37 (6H, m), 1.24 (38H, br s), 1.04-1 00 (9H, 3 s), 0.84 (3H, tr, J 
= 6 .56 Hz) FAB-MS [M+H ] calcd for Cj^NjOgPSi 943 4, found 943.0. 

l-Myristoyl-2-(13-carboxytridecanoyl)-s/i-3-phosphatidylglycerol-2 , -3'- 
isopropylidene (la, PG J ^°° hI/p ) 

30 

...... , v „.va.L„ i c suspension was <;Tirrea < (_ : - m ; r , ■ 

before pipetting 10 62 ml of NaOAc buffer (composed of 100 mM NaOAc and 50 mM 
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CaCk pH 6.5) and 1 6 ml (12 10 mmol) of (R)-isopropylidene glycerol (IPG) After 
addition of IPG, the two-phase system immediately emulsified PLD (100 \xl 1.0 
unit/iil in NaOAc buffer) was added and the reaction monitored by TLC 
(CHC1 3 CH 3 OH NH 3 H 2 0 = 65 :25 :4), after 2 h > 60% of the PC-COOH was 
5 converted to pG 2taHCOOKP . Additional PLD (40 ^1) was added and the reaction 

incubated for another 5 hours (90 % completion based on TLC). The reaction mixture 
was acidified using KHS0 4 (1 M, 1.0 ml), then concentrated by rotoevaporation and 
the resulting residue was purified by flash chromatography using CHC1 3 /CH 3 0H/H 2 0 = 
65:25:3. Chromatography fractions containing product were pooled, concentrated by 

10 rotoevaporated and vacuum dried; the yield after column chromatography was 1.74 g 
of a white solid (83.2% yield) TLC in CHC1 3 /CH 3 0H/H 2 0 (b5:25:4) showed a pure 
product (Rf = 0.40). FTIR (thin film, cm' 1 ) 2982.51, 2923.36, 2852.97, 1741.51, 
1710.84, 1556.64, 1462.17, 1415.60, 1379.37, 1370.00, 1244.68, 1219.69, 1160.35, 
1 106.00, 1069. 15. J H NMR (500 MHz, CDC1 3 , 6 ppm): 5.24 (1H, br, m), 4.37 (1H, 

15 br, m), 4.29 (1H, br, m), 4.17 (2H, br, m), 4.04 (3H, br, m), 3.90 (1H, br, m), 3.80 
(1H), 2.29 (6H, m), 1.57 (6H, m), 1.42 and 1.35 (3H, s, -CCH 3 , one of the two 
conformers), 1.40 and 1.33 (3H, s, -CCH 3 , one of the two conformers), 1.26 (36H, m), 
0.86 (3H, t, J = 7 00 Hz) FAB-MS [M+H + ] calcd for C 37 H 69 0 12 P 737.9, found 737.5. 

20 l-myristoyI-2-(13-carboxytridecanoyI)-s/i-3-phosphatidyl-p-nitrophcnyIcthanol 
(lb, PA 2 "- COOH/p ) 

PA 2u ^°° HP was prepared similar to PG 2u ^°° HP except that p- 
nitrophenethyl alcohol was used instead of IPG for the PLD transphosphatidylation 

25 reaction and the reaction was scale up by 15 fold Optimum yields required PLD (1.0 
unit/nl) to be sequentially added in three increments 100 (al, 2 h incubation; 37 jal, 2 h 
incubation, then 50 nl-PLD, incubation for 12 hours). The product was both phospray 
positive and UV positive The reaction mixture was acidified with KHS0 4 , (1 .042 ml 
of 1M), concentrated by rotoevaporation and the product purified by flash 

30 chromatography usine silir^ t?M TV -< 1 * • 

ontainmg product v. ere poo^ci. ^uiiL.eniraiea o> rotoevaporated and vacuum dried, 
the yield after chromatography was 1 12 g of light yellow solid (75. 1% yield) TLC in 
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CHC1 3 /CH 3 0H/H 2 0 (65:25:4) showed a pure product (Rf = 0 42). FTER: (thin film, 
cm' 1 ) 2923.70, 2852.81, 1740.44, 1705.13, 1603.81, 1550.54, 1522.1 1, 1465.95, 
1345 88, 1220 07, 1108.70, 1070 67, 1034.59. FAB-MS [M+H*] calcd for 
C 39 H^N0 12 P 772.9, found 777.5 

5 

N-/err-Butyloxycarbonyl-l-myristoyl-2-(13-carboxytridecanoyI)-5/i-3- 
phosphatidylethanolamine (lc, PE 2 ^ 00 " 7 *) 

PE 2u^oohp was preparec j similar to PG 2u ^°° HP except that a mixed 

10 solvent system, ethyl acetate/chloroform (4.6:1) was used Chloroform was alcohol 
free Based on scanning densitometry of TLC plates developed in 
CHC1 3 :CH,0H:NH 3 H 2 0 (6V25:4), -88 0% of PC-COOK starting material was 
converted to the corresponding PE carboxyl analog. The reaction mixture was 
concentrated by rotoevaporation and the intermediate was not purified prior to 

15 protecting the primary PE amine with Boc. To the concentrated residue, Na 2 C0 3 (20 
ml, 1 .0 M) and dioxane (40 ml) were added. The reaction mixture was cooled to 0 °C 
before adding (Boc) 2 0 (1 98 g, 9 07 mmol). The reaction was stirred (0 °C, 1 h), r.t. 
for 3 hours followed by acidification by KHS0 4 (40 ml, 1 .0 M) Addition of 
chloroform (40 ml) caused the milky emulsions reaction mixture to form two layers. 

20 The organic layer was collected and combined with the chloroform washings (2 x 40 
ml) of the aqueous layer The chloroform extractions were concentrated by 
rotoevaporation before purification of the product using silica gel flash 
chromatography. The solvent system was CHCl 3 /CH 3 OH (9: 1, 300 ml) followed by 
CHC1 3 /CH 3 0H/H 2 0 (65:25:3, 800 ml) Chromatography fractions containing product 

25 were pooled, concentrated by rotoevaporated and vacuum dried, the yield after 
chromatography was 0 96 g of white solid (88 .8 % yield) TLC in 
CHC1 3 CH 3 OH NH 3 H 2 0 (65:25:4) showed a pure product (Rf = 0.44) FTER (thin 
film, cm" 1 ) 3384.65, 2923 00, 2852.32, 1742.15, 1715.33, 1690.96, 1553.86, 1517.93, 
1464 04, 1365.55, 1244.64, 1221.96, 1171 99, 1108 70, 1071.16 ! H NMR (500 

30 MHz, CDCL 6 ppmV 5 20(1H br ^ 4^nu ^ t 1 ^ 

. iwi. >. . - . rv- ; ^ ; : h . : ~ ; \w > t At* -Mb [M*Na J eaicU 

for C 3g H 71 N0 12 P 789 0, found 788 5 
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N-/err-But>loxycarbonyl-)-aIlylester-l-myristoyI-2-(13-carboxytridecanoyl)-5/i-3- 
phosphatidylserine (Id, ps 2 ^°° HT ) 

PC-COOH (7.0 g, 9.90 mmol) was dissolved in anhydrous chloroform 
5 (200 ml) using a 500 ml r.b. flask. Serine allyl ester (7.2g, 49.50 mmol) was added to 
the reaction mixture followed by 33 .4 ml of buffer (composed of 100 mM NaOAc, 
50mM CaCl 2 at pH 6.5). PLD (500 n!, 1.0 unit/jil in NaOAc buffer) was added and 
the reaction was stirred at 30 °C under N 2 for 8 hours. The chemical intermediate, a 
carboxyl PS analog from the PLD reaction of PC-COOH, was concentrated via 

10 rotoevaporation and not purified prior to the protecting the serine amine with Boc 
The primary serine amine was protected with Boc as described above in the 
preparation of PE 2u ^ OOHP After acidification of the protection reaction mixture using 
KHS0 4 (160 ml, 10 M), the product was extracted using diethyl ether (4 x 100 ml) 
The diethyl ether extracts were pooled and the solvent rotoevaporated to obtain a 

1 5 yellow oil which was purified by silica gel flash chromatography using 

CH 2 C1 2 /CH 3 0H/H 2 0 (65:25 :4) as the mobile phase. The amount of final product 
(white solid) obtained after column chromatography resulted in an 87% overall yield 
(7.32 g) TLC in CH 2 Cl 2 /CH 3 OH (1:1) showed a pure product (Rf = 0 81) FTIR: (thin 
film, cm" 1 ) 3364.26, 2925.48, 2855.87, 1741.78, 1713.99, 1510.08, 1461.08, 1372.02, 

20 1185.02, 1080.51. ] H NMR (300 MHz, CDC1 3 , ft ppm): 5.90 (1H, m), 5.34-5.20 (4H, 
m), 4 .62-4 .00 (8H, m), 2.27 (6H, m), 1.57 (6H, m), 141 (9H, s), 1.25 (36H, br s), 
0.87 (3H, tr, J = 6 65 Hz). MALDI-MS [M-H]" calcd for C 42 H 75 N0 14 P 848.1, found 
848.0. 

25 " trr IAM.PG c,0/CJ (7) 

PG 2u-coohp (i a ) (i.o g, 1 35 mmol) was suspended in anhydrous CHC1 3 
(50 ml) using a r.b flask (100 ml) and CDI (0.380 g, 2.35 mmol) was added. The la 
Hgand is not very soluble in CHC1 3 , but, the CDI-activated pg 2 ^ 00 ^ ligand is very 
soluble After overnight stirring at r t , both TLC CH 2 CyCH 3 0H/H 2 0 (65:25:4) and 

30 IR spectroscopy indicated that CDT ~ — " - J 1 

tne suspension \kx>* mixec using an oruitai Miaker { i 5u rpm, J 4 hours) The 
""lAM PG T ' silica was then filtered through an F glass-sintered funnel, washed with 
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chloroform (2 x 20 ml), methanol (3 x 20 ml) and acetone (10 ml) followed by drying 
in a vacuum oven at 40 °C overnight The coupling was repeated using reaction 
conditions corresponding to 0 25 equivalents of the pg 2u>jCC>oht> ligand. C10 and C3 
endcappings were performed as described previously to obtain c * tCT IAM.PG pcl0C3 . 
5 Deprotection of the isopropylidene protecting group on the e " tff IAM.PG p ciaC3 surface 
was needed to generate an immobilized artificial membrane containing PG. Briefly, 
" to lAM.PG p/CI0/C3 (-10.9 g) was suspended anhydrous CU 2 C\ 2 (24 ml). Then 
anhydrous TFA (8 ml) was added to the suspension. After 2 h of shaking the silica 
was filtered, washed with CH 2 C1 2 (2 x 20 ml), CH 3 OH (3 x 20 ml), and acetone (20 

10 ml), followed by vacuum drying (40 °C overnight). The deprotection step was 

repeated with 50% TF A/50% CH 2 CK (V/V, 24 ml) for 4.0 hours. Then the silica was 
immersed in 50% 1 M Na 2 CO 3 /50% CH 3 OH (80 ml, 6 min) under constant stirring, 
followed by washing with water (600 ml), CH 3 OH (2 x 40 ml) and acetone (2 x 40 ml). 
The surface was finally vacuum-dried at 40 °C overnight. FTIR microscopy of surface 

1 5 (7) indicated that the IPG group was quantitatively removed because the characteristic 
IPG peak at 1370.0 cm" 1 disappeared completely after the deprotection. The bonded 
ligand density[Ong, 1994 #14] (la) on ertCT IAM.PG p was 134.67 nmol/g-silica based on 
elemental analysis of carbon content before (C, 2.13%) and after (C, 7.46%) ligand la 
bonding to silica surfaces. 

20 

«"TAM.PA C10/C3 (8) 

The bonding procedures for immobilizing pa 2 ^ 00 ^ (lb), C10 and C3 
endcappings to prepare e * tCT lAM.PA pcl0C3 were the same as those used to prepare the 
wtCT IAM PG pcl0 C3 Deprotection of p-nitrophenylethyl (NPE) group on the 
25 nAM.PA pc,on surface was performed as described Briefly, e * tCT lAM.PA pcl0C3 

(-15.0 g) was suspended in anhydrous CH 2 C1 2 (50 ml), then anhydrous DBU (20 ml) 
was added After 2 h of shaking, the silica was filtered, washed with saline (1.0 M 
NaCl, 500 ml), H 2 0 (500 ml), and acetone (50 ml), followed by vacuum drying (40 °C, 
overnight) Based on the integrated FR band of VPF 1 • 7 " 1 ^' c 

. . '. i 'pnenviet fi\ ; croup \ m\ . >^ .■ deprotecttop using DC-nzenc ic^uiteo in - " o 
After deprotection. the silica was filtered, washed with CHCl 3 (200 ml). CH 3 OH (100 
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ml), saline (10 M NaCl, 500 ml), H 2 0 (500 ml) and acetone (100 ml) followed by 
vacuum drying (40 °C overnight) The bonded ligand density (lb) on c * tCT IAM.PA p was 
96 50 lamol/'g-silica based on elemental analysis of carbon contents before (C, 2.13%) 
and after (C, 6.34%) ligand lb bonding to silica surfaces. 

5 

- U TAM.PE C10/C3 (9) 

The bonding procedures for immobilizing pe 2 ^ 00 ^ (lb), C10 and C3 
endcappings to prepare ortCT I AM . PE p/c 1 0/03 were the same as those used to prepare the 
"■"I AM . P G p c 1 o c3 . Boc Deprotection from the ""^T AM . PE p/c 1 0 C3 surface was 

10 performed as described Briefly, c * tCT I AM . PE P/C 1 0/C3 (-8.9 g) was suspended in 
anhydrous CH 2 C1 2 (20 ml). Then anhydrous TFA (20 ml) was added to the 
suspension. After 2h of shaking, the silica was filtered, washed with CH 2 C1 2 (2 x 20 
ml), CH 3 OH (3 x 20 ml), and acetone (20 ml), followed by vacuum drying (40 °C, 
overnight). The first deprotection was - 80 % complete and a second deprotection 

15 step in 50% TF A/50% CH 2 C1 2 for 4.0 hours resulted in complete deprotection, the 
Boc IR band 1369.0 cm' 1 completely disappeared. The bonded ligand density (lc) on 
C * CT IAM.PE P was 109.33 |jmol/g-siIica based on elemental analysis of carbon contents 
before (C, 2.13%) and after (C, 6.64%) ligand lc bonding to silica surfaces. 

20 " U TAM.PS CI0/C3 (10) 

The bonding procedures for preparing CTtCT IAM.PS p C3 C3 were the same 
as those used to prepare the ^AM.PG^ 01003 except that double C3 endcapping was 
used Deprotection of surface allyl ester group on c " tCT IAM.PS P C3 C3 was based on the 
solution deprotection conditions Briefly, e * tCT lAM PS P C3C3 (22 44 g) was suspended in 

25 anhydrous CH 2 C\ 2 (50 ml), nitrogen purged (5 min), sonicated (10 sec) and nitrogen 

purged (1 min) before adding the deprotection reagents. PdCl 2 (PPh 3 ) 3 (96.07 mg, 0. 14 
mmol), acetic acid (600.4 jil, d=1.049, 10.50 mmol) and Bu 3 SnH (4.53 ml, d- 1.082, 
16 84 mmol) were sequentially added to ~ tCT IAM PS P ' C3 C3 suspension After 2h of 
shaking, Bu,SnH (1 0 ml. 3 72 mmol) was added and *he v>mf>~^ — • - • v ' 



:v:\anc 1 ; rr, . > t 1 * : ;e:ai^' i ; . mi : ana acetone < 4* - nu ), ioiiowed D\ 
vacuum drying (40 C C, overnight) The deprotection procedure of the allyl ester wa> 
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repeated prior to Boc deprotection For Boc deprotection, the conditions described 
above to prepare the C * ICT 1AM.PE C10C3 surface were used except that the deprotection 
reaction was performed once. Based on FTIR microscopy, quantitative deprotection 
occurred The bonded ligand density on ertcr IAM.PS p was 104.86 ^mol/g-silica based 
5 on elemental analysis of carbon contents before (C, 2. 13%) and after (C, 6.99%) ligand 
Id bonding to silica surfaces 

IAM.SM C10C3 (11) 

The bonding procedures for immobilizing SM^ 00 * 1 * (2), CIO and C3 
10 endcappings to prepare IAM SM P C10 C3 were the same as those used to prepare the 

wlcr IAM.PG p C10C3 except the ligand coupling to silica propylamine was performed only 
once Deprotection of TBDPS protecting group on the SM headgroup was done as 
described. Briefly, IAM.SM PC!0C3 (-2.1 g) was suspended in THF (20 ml) before 
TBAF (20 ml, 1 .0 M solution in THF) was added to the suspension. After 2 h of 
15 shaking, the silica was filtered, washed with THF (2 x 20 ml), CHC1 3 (3 x 20 ml), 
CH 3 OH (3 x 20 ml), and acetone (20 ml), followed by vacuum drying (40 °C, 
overnight). FTIR microscopy indicated that deprotected was quantitative, because the 
TBDPS IR peak at 3074.23 cm' 1 and 3056.10 cm" 1 completely disappeared. The 
bonded ligand density (2) on IAM SM P was 90.53 ^mol/g-silica based on elemental 
20 analysis of carbon content before (C, 2.04%) and after (C, 7.35%) ligand 2 bonding to 
silica surfaces. 

Synthesis and Immobilization of Ceramide Analogs on the Silica Particles 

The skin is the largest organ of the human body, a primary area 
25 contacted with the environment, and a route by which many chemical substances enter 
the body Research has demonstrated that drug delivery through the skin is feasible for 
many simple potent drug molecules (less than 1000 Da) via transdermal drug delivery 
systems Seven drugs marketed in United States including clonidine, estradiol, 
fentanvl nicotine nitrnelvrenn s.- - ^ M • 1 ' * 

_ , irk-enirations o! a a rug can f-e d^me\eu \utnout the high peaK levels associated with 
oral therap\ The avoidance of high peak levels may help minimize the side effects of 
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certain drugs Meanwhile dermal exposure to toxic substances represents a major 
occupational hazard, successful anticipation of potential risk could significantly reduce 
the incidence of chronic health and environmental problems. The interest in the skin 
has encouraged research in the field and has led to a better understanding of skin 
5 biological structure The studies have showed that the major permeability barrier of the 
human skin is provided by its outer layer - the stratum corneum (SC). The stratum 
corneum consists of dead cells surrounded by an extracellular matrix containing lipid 
lamellae The major lipid components of this extracellular matrix is ceramides which 
comprise 50% of the total lipid and consist six structurally heterogeneous ceramides 

10 Among them, ceramide 2 comprises 40% of the total ceramides. Ceramide consists of 
mainly 24- through 28- carbon fatty acids amide-linked to sphingosine and 
dihydrosphingosine bases. 

Studies of model systems capable of predicting plasma levels following 
topical administration are continuously being pursued. These investigations aim at 

15 facilitating the rational selection of transdermal drug candidate which can penetrate 
skin barrier or predicting the potential risk of dermal exposure to toxic chemicals. By 
immobilizing skin lipids such as ceramide 2 on a chromatographic support there is 
provided a new model system to evaluate chemical skin interaction. The chemical skin 
interactions can then be chromatographically measured via HPLC. This system would 

20 be much better than octanol/water partition system and more feasible than laborable 

and controlled measurement of skin permeability in vivo. The ceramide ligand (N-[13- 
carboxyltridecanoyl]-D-eryf/?ro sphingosine 2) used to prepare ceramide silica surface 
contains basic ceramide 2 functional moiety, and a free o> -carboxyl group. The free a) 
-carboxyl group functions as a tag to the silica propyl amine particles through 1,1'- 

25 carbonyldiimidazole (CDI) activation Anhydrides are used to endcap excess amine 
groups on the silica propyl amine particles See Scheme 5. 
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A representative synthesis of ceramide stationary surfaces is outlined in 
Scheme 6 The treatment of N-Boc-L-serine with methyl iodide gave N-Boc-L-serine 
methyl ester 3 which was reacted with 2,2-dimethoxypropane under the catalysis of p- 
toluenesulfonic acid monohydrate to produce the acetal protected compound 4, which 
5 was reduced to the aldehyde 5 under mild reduction (using diisobutylaluminum hydride 
(DIBALH) as reagent and running the reaction at -78°C). Meanwhile, 1-pentadecyne 
6 was synthesized from 1-bromotridecane and lithium acetylide ethylenediamine 
complex in dimethyl sulfoxide (DMSO) solution Then oxazolidine aldehyde 5 was 
treated with lithium acetylide (from 1-pentadecyne 6 and n-butyl lithium) at -23 °C to 

10 produce the basic structure of sphingosine 7 in 84% yield after flash chromatography. 
High erv7/im-se!ectivity (with diastereoseieciivity (ds) - 89%) was observed in this 
reaction. Treatment of compound 7 with mild acid Amberlyst 15 in methanol at room 
temperature resulted in selective cleavage of the acetal protecting group, leading to the 
1,3-diol 8 which was converted to N-Boc-sphingosine 9 by selective reduction of the 

15 triple carbon bond with Red-Al in ether. This reduction condition specifically 

converted the triple bond to the trans double bond. Cleavage of the carbamate moiety 
with 1.0 N HC1 in dioxane at 100°C resulted in D-ery^ro-sphingosine 10 which was 
reacted with tert-butyldiphenylsily chloride (TBDPSC1) in order to protect the two 
hydroxyl groups The TBDPS group is considerably more stable (about 100 times) 

20 than the tert-butyldimethylsilyl (TBDMS) group toward acidic hydrolysis and therefore 
is appropriate for the immobilization reaction. Subsequent reaction of TBDPS- 
protected sphingosine 11 with dodecanedicarboxylic acid anhydride 12 (a cyclic 
anhydride, synthesized from diacid and dicyclohexylcarbodiimide (DCC)) in 
chloroform afforded N-[13-carboxyItridecanoyl]-D-er>7/ira-sphingosine 2 under 

25 catalysis of 4-dimethylaminopyridine (DMAP) in 87% yield after flash 

chromatography The o>-carboxyl group of ceramide analog 2 was activated with CDI 
in chloroform Following the activation, the ceramide-imidazolides were bonded to 
silica propylamine (SPA) particles with 24 hours shaking After workup, ceramide- 
based silica sta'ion^rv n^o- ■> 1 * '■ ' 

-moved wjtr, ;ctraDui> idminonium liuonde(TBAF) in THF The reactions were 
monitored by a Nicolet Magna 550 FT-IR spectrometer equipped with a Spectratech 
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IR-Plan I microscope. After deprotection, a monolayer of ceramide lipid membrane is 
generated bearing identical interfacial groups to endogenous ceramide lipid. The 
bonding density was calculated from elemental analysis. Finally these silica particles 
were packed into the stainless-steel columns and used for high performance liquid 
5 chromatography Our preliminary data clearly showed that ceramide based stationary 
phases are suitable for the prediction of skin permeability constants of diverse 
chemicals Solute capacity factors (k') measured on the columns correlated well with 
skin permeability coefficients (kp) or the percentage of percutaneous absorption (% 
Abs) measured through the excised human skin or the human subjects. 

10 

Correlation of Drug Bi nding Affinities with Oral Absorption 

Table 6 shows eight compounds eliciting low oral absorption that were 
used as a training set to calculate a tow MAF u representative of low drug absorption. 
The small average D 2 value, i.e., 3.5 ± 1.24 relative to k)W MAF fi (last column in Table 

15 6) indicates that the group of compounds in the training set are statistically well 

defined. D 2 values based on lovs MAF M for the 13 compounds exhibiting complete oral 
absorption are very high (Table 7). The large average D 2 value, i.e., 96.26 ± 96.01 
(last column in Table 7) indicates that compounds which are highly absorbed have 
statistically different membrane binding properties compared to compounds that exhibit 

20 low oral absorption. 
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Table 6 Training Set Representing Low Oral Absorption 



Compound 


Oral Absorption (%) 


D : ** 


nalbuphine 


16 


2.10 


ketamine 


20 


3.71 


scopolamine 


27 


2.13 


norfloxacin 


30 


4 96 


sulpiride 


36 


4.95 


acebutolol 


37 


2 87 


ranitidie 


50 


4.67 


ptr*nr*ciHg 
j -~ 


cn 


2.58 






Average - 3 .5 ± 1.24 



♦These eight compounds were used to determine the kw MAF M vector. 
**D 2 values are calculated relative to the ^MAF* 1 vector 

15 
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Since the test compounds in Table 7 exhibited distinct D 2 values 
compared to the km MAF M vector, the compounds listed in Table 7 were used as 
potential training set to generate a ^MAF^ 1 vector 

Table 7. Calculation of D 2 Relative to km MAF M for 



urn puunu 


Oral Ahcnrnfinn (%%\ 


D 2 


uumuinc 




22.4 


cionoraic 




a, J "™r 


dlUI CI1U1U1 


\j 


27.8 


m^t rnntH q 7nlp 

llltll VJl Li L1CL^.VJ J C 


98 


81 8 


i n f\ r\ m p t h a r i n 


98 


65 3 


clonazepam 


98 


151 


desmethyldiazepam 


99 


328 


phenytoin 


100 


7.17 


phenobarbital 


100 


17.3 


diazepam 


100 


138 


corticosterone 


100 


17.2 


chlordiazepoxide 


100 


103 


caffeine 


100 


58.4 




Average = 96.26 ± 96.01 
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Table 8 shows D 2 values relative to lu8h MAF M for each compound in the 
high drug absorption training set D 2 values based on ^MAF* 1 for the 8 compounds 
exhibiting low oral absorption are high (Table 9). 

5 Table 8. Training Set Representing High Oral Absorption 



Compound 41 


Oral Absorption (%) 


J mm 

D 


chlordiazepoxide 


100 


1.26 


diazepam 


100 


1.67 


clofibrate 


95 


1 85 


clonazepam 


98 


1 89 


alprenolol 


96 


2 48 


phenobarbital 


100 


2 99 


corticosterone 


100 


3 05 


phenytoin 


100 


3 16 


chlonidine 


95 


4 87 


desmethyldiazepam 


99 


5.57 


indomethacin 


98 


5.75 


caffeine 


100 


6.22 


metronidazole 


98 


7.15 




Average = 3 .68 ± 1.97 



*These compounds were used to determine the ^MAF^ vector 
**D 2 values are the calculated relative to the * ugh MAF M vector. 



25 
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Table 9 Calculation of D : Relative to h * h MAF M for Compounds 



1 


nat Have Low Oral Absorption 


€^ n in r*n ii n rl * 


Oral Absorption (%) 


D 2 ** 


nalhnnhinf* 

J 1AJ ULi Lf I Lll It 


16 


2.20 


NCI dl 1 LI 1 1 fc 


20 


7.87 


cp orvol 3 rn i n f* 


27 


2.53 


norfloxacin 


30 


39 4 


sulpiride 


36 


3.94 


acebutolol 


37 


2.05 


ranitidine 


50 


42.09 


— : j - 

VlULiUSlUC 


50 


3.36 




Average - 12.9 ± 17.3 



15 Data Acquisition Using Liquid Chromatograph/Mass Spectrometer (LCMS) - 
General Experimental Protocol 

Mass spectrometers (MS) can be used as detectors for identification of 
compound eluting from a liquid chromatograph (LC) instrument The only 

20 requirement for detection is that the compound can be ionized without composition. 
Thus, one preferred method for determining membrane binding constants is to use an 
MS as a detector for compounds eluting from LC columns. The general advantages of 
using MS detection are (1) compounds without chromophores can easily be detected, 
(2) multiple compounds can be uniquely identified from a single injection, and (3) small 

25 amounts of compound can be used for determination of binding constants 

chromatographically measured on immobilized artificial membranes (IAMs) or other 
solid phase chromatographic substrates. These three advantages allow simultaneous 
injection of very small amounts of hundreds of compounds with detection of essentially 
all compounds eluting from the column Such procedure allows for direct comparison 

30 of binding affinities of a large number of test compounds with multiple compound* 

. ,,. .. .j; ■. . .an;.;: \ an,: wiTiannir, n 1 experimental 
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parameters. An example of data acquisition using LC with MS detection is given 
below 



Test Mixture 

5 A mixture containing 0.100 jamoles of 155 compounds was prepared. 

At least a portion of the compounds are compounds of known clinical significance, 
they serve as internal standards against which the binding affinities of the test 
compounds are measured. The mixture was dissolved in 6.0 mL 15% acetonitrile and 
85% 30mM ammonium acetate, pH - 7.4. 20 (iL of this test mixture was loaded onto 
10 a 3 x 0.46 cm I AM column This corresponds to loading ~ 16 pg of the test mixture 
and loading - 320 pmoles of each compound onto the column 



LC Conditions 

A 3 x 0 46 cm IAM column was used. The mobile phase was 15% 
15 acetonitrile and 85% 30mM ammonium acetate, pH - 7.4. The flow rate was 1.0 mL 
per minute and no flow splitter was used. The run time was 2 hours followed by a 10 
minute 100% acetonitrile wash. 



MS Conditions 

20 A Bruker Esquire-LC was used Each run was done twice. In the first 

run, the auto MS/MS was not used In the second run, the auto MS/MS was used 
The auto MS/MS allows for the identification of compounds which have the same 
molecular weight 

Typical MS instrument have the capacity to use flow rates of 2 mL/min 
25 and at least one instrument has a mechanism for washing salts from the inlet of the ion 
source Typically, mobile phases cannot have salts because evaporation of the mobile 
phase in the ion source results in excess buildup of salts that restrict flow of ions into 
the ion source Finnigan recently developed an LC/MS instrument in which salts are 
washed away during the chromatography process Affinitv data are collerteH ^ 

.v.^l^ . a nipuic; a a: a accessinie ,;om rr.ur.i cat icr. the udid 
storage device, and programmed with a vector analysis algorithm, to provide a user 
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readable output in a preprogrammed format for facilitating comparison of the binding 
affinities of test compounds with those of control compounds or with mean vector 
quantities for one or more training sets 



